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Using  both  the  plane-wave  and  the  Korringa-Kohn-Rostoker  method,  the  photon  band  structure 
is  calculated  within  the  scalar  wave  approximation  for  a  fee  lattice  of  dielectric  microspheres  em¬ 
bedded  within  a  uniform  host  medium.  The  plane-wave  method  is  found  to  converge  fairly  rapidly. 
The  optimal  volume-filling  fraction  of  spheres  for  the  creation  of  a  gap  in  the  photon  density  of 
states  is  found  to  vary  substantially  with  the  relative  dielectric  constant  r,  and  this  gap  persists  for  r 
as  small  as  3. 


There  has  been  growing  interest  in  studying  photon 
band  structures  in  three-dimensional  periodic  dielectric 
media  that  exhibit  a  gap  in  the  photon  density  of 
states.1-3  There  are  important  motivations  for  such 
studies.  (1)  Since  electromagnetic  modes  are  totally  ab¬ 
sent  within  the  band  gap,  spontaneous  emission  is  there¬ 
fore  strongly  suppressed.  This  ability  to  inhibit  spontane¬ 
ous  emission  has  some  extremely  important  consequences 
since  spontaneous  emission  plays  a  fundamental  role  in 
limiting  the  performance  of  many  optical  and  electronic 
devices  such  as  semiconductor  lasers,  heterojunction  bi¬ 
polar  transistors,  and  solar  cells.1,4  In  atomic  physics, 
this  phenomenon  has  in  fact  been  demonstrates.5-7  (2)  It 
has  also  been  pointed  out  that  the  basic  properties  of 
many  atomic,  molecular,  and  excitonic  systems  can  be 
profoundly  modified  in  a  volume  of  space,  where  the 
most  important  electromagnetic  processes  are  totally  ab¬ 
sent.8  (3)  There  is  also  a  proposal  to  study  mobility  edges 
and  Anderson  localization  of  photons  within  a  related 
pseudogap  in  the  presence  of  some  randomness.2 

The  first  calculation  of  photonic  bands  exhibiting  the 
presence  of  a  common  gap  in  the  photon  density  of  states 
was  carried  out  recently2  based  on  the  scalar  wave  equa¬ 
tion  and  using  the  Korringa-Kohn-Rostoker9,10  (KKR) 
method.  The  structure  consists  of  a  fee  lattice  of  micro¬ 
spheres  embedded  in  a  host  medium.  The  main  results  of 
that  study  were  (1)  the  gap  in  the  photon  density  of  states 
persists  down  to  a  relative  dielectric  constant  r  of  about 
7.8,  and  (2)  the  optimal  volume  filling  fraction  /,  for  the 
creation  of  this  gap,  is  about  11%.  In  this  work,  we  re¬ 
port  results  that  are  based  on  the  plane-wave  method. 
Results  based  on  the  KKR  method  are  also  presented. 
We  find  that  our  plane-wave  calculation  converges  fairly 
rapidly.  This  is  an  important  finding  since  the  plane- 
wave  method  is  comparatively  much  simpler  than  the 
KKR  method,  the  computer  program  is  much  easier  to 
write,  and  runs  substantially  much  faster.  Moreover,  un¬ 
like  the  KKR  method,  which  is  limited  to  “muffin-tin”- 
like  modulation  of  the  “potential,”  the  plane-wave 
method  can  readily  handle  all  sorts  of  modulations.  This 
capability  is  especially  important  since  we  believe  that, 
besides  requiring  a  high  dielectric  contrast,  the  shape  of 
the  dielectric  structure  is  also  a  very  crucial  factor  in 


determining  the  photon  band  structure  and  the  existence 
of  a  gap  in  the  density  of  states. 

Although  some  of  our  results  are  somewhat  similar  to 
those  reported  earlier,  there  are  two  significant  discrepan¬ 
cies.  First,  our  results  show  that  the  gap  persists  done  to 
an  r  value  of  about  3.  Second,  we  find  that  the  optimal 
value  of  /  depends  significantly  on  the  value  of  r,  and  in 
general  decreases  with  increasing  r. 

Following  the  previous  work,2  we  start  with  the  scalar 
wave  equation 

~y2\p+Vip—klrij ,  (i) 


where  the  “potential”  V  is  given  by 


V=kt 


(2) 


with  e  —  ea  inside  the  spheres,  e~eb  inside  the  host,  and 
kb  =  y/ nebco/c.  Although  the  equation  is  of  scalar  form, 
and  therefore  neglects  the  vector  nature  of  the  photon,  it 
retains  an  important  characteristic  in  that  the  “potential” 
is  proportional  to  <u2,  and  thus  vanishes  in  the  long- 
wavelength  limit.  This  has  some  very  important  physical 
consequences  in  the  present  problem  as  well  as  in  the 
photon  localization  problem.11 

The  plane-wave  method  is  very  straightforward.  We 
work  in  Fourier  space,  where  the  Fourier  coefficients 

ck  =  j*dre-,kr0(r)  ,  (3) 

and 

Vk  =  ~nf  dre~iK  rV(r)  ,  (4) 

fi  cell 

with  K  as  a  reciprocal-lattice  vector  of  the  fee  lattice,  and 
fi  the  volume  of  the  unit  cell.  The  scalar  wave  equation 
can  then  be  expressed  in  the  form 

[*ft2Hk  — K!2]ck_K-*22UK_K,ck_K.=0  ,  (5) 

K 

where  we  have  defined  U—V /kb.  This  eigenvalue  equa¬ 
tion  can  be  recast  into  the  standard  form  if  we  divide  it 
by  kb\k  —  K|,  let  dk_K  =  |k  —  K|ck_K,  and  treat  k^2  as 
the  eigenvalue.  The  result  is 
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1-t/p 

|k-K|2 


_1_ 

kl 


<^k-K 


I 


Uk-k 

Ik-K'| 


d k- 


K' 


=  0  ,  (6) 


J  where  the  prime  over  the  summation  sipn  signifies  that 
the  K'  =  K  term  is  to  be  omitted.  The  photon  band  struc- 
j  ture  is  then  obtained  by  solving  this  equation  for  the  ei- 
I  genvalues  kf2  for  each  value  of  k. 
j  Before  we  give  the  results  of  our  calculations  that  are 
'  specific  to  the  present  problem,  we  want  to  make  a  few 
]  general  remarks  on  the  band  structure  based  on  the  above 
;  plane-wave  band  equation.  First,  it  is  clear  that  the  equa- 
i  tion  has  exactly  the  same  form  as  the  corresponding 
j  equation  for  electrons,  if  we  identify  (1  —  C/0)/|k  —  K|2 
!  with  the  free-particle  energy,  kb2  as  the  exact  energy  ei- 
5  genvalue,  and  I/K_K-/(|k  — K||k— K'|)  as  the  corre- 
;  sponding  Fourier  element  of  the  electron  potential. 

'  Thus,  in  the  empty-lattice  limit,  i.e.,  V— *0,  the  band 
;  structure  can  be  obtained  from  most  solid-state  text- 
■  books.  It  is  clear  that  in  this  limit,  most  of  the  levels  are 
;  highly  degenerate,  especially  at  high-symmetry  points, 
1  and  for  k  varying  from  the  T  point  to  the  edge  of  the 
i  Brillouin  zone,  the  dispersion  curves  are  straight  lines 
j  given  by  kb  =  k. 

j  There  are  two  important  effects  when  the  “potential”  is 
!  turned  on.  First,  depending  on  the  symmetry  of  V,  some 
I  of  these  degeneracies  are  lifted.  Second,  the  levels  shift 
]  up  away  from  zero  if  U0  is  positive,  and  shift  down  to- 
I  ward  zero  otherwise.  As  a  result,  the  dispersion  curves 
originating  from  the  T  point  are  linear  only  near  the  T 
point,  where  k  is  small  compared  with  the  magnitude  of 
the  smallest  K  of  the  lattice,  and  the  slope  of  the  straight 
portion  is  no  longer  unity,  but  should  be  given  by 
(efc/eefr),  where  eeff  is  the  effective  dielectric  constant  of 
the  entire  medium  in  the  long-wavelength  limit.  The  ex¬ 
act  expression  for  eefr  can  be  easily  derived  from  the 
plane-wave  band  equation  because  at  low  frequencies  and 
long  wavelengths,  where  k£  and  k  are  both  small  com¬ 
pared  with  the  smallest  nonzero  reciprocal-lattice  vector, 


r-1125,  f-0.373 


FIG.  1.  The  solid  curves  give  the  band  structure  calculated 
with  the  plane-wave  method.  The  volume-filling  fraction  of 
spheres  is  0.375  and  the  relative  dielectric  constant  is  12.25. 
The  solid  circles  are  our  results  calculated  with  the  KKR 
method. 


FIG.  2.  The  optimal  volume-filling  fraction  of  spheres  for  the 
creation  of  a  gap  in  the  photon  density  of  states  is  plotted  as  a 
function  of  the  relative  dielectric  constant. 


only  the  coefficient  dk_K  with  K=0  is  non-negligible. 
Thus  one  can  easily  show  that 

ceir=(  1  —  U(>)eb  •  (7) 

Now  we  are  ready  to  present  our  results  for  the  present 
model.  For  the  case  of  dielectric  spheres  considered  here, 
we  have 


UK=3f(l-r)g(Ka)  ,  (8) 

where  the  function  g  ( x )  =  ( sin*  —  x  cosx )  /x 3,  a  is  the  ra¬ 
dius  of  the  sphere,  and  U0=f{  1  —  r).  Using  these  results 
and  the  plane-wave  equation,  we  have  calculated  the  pho¬ 
ton  band  structures  for  various  values  of  /  and  r.  We  find 
that  the  results  for  the  lowest-lying  bands  converge  fairly 
rapidly.  To  within  an  accuracy  of  0.1%,  we  find  that  70 
K  points  are  sufficient  for  convergence  even  for  r  as  large 
as  12.  For  r  around  3,  only  30  points  are  needed.  This 
rapid  convergence  is  reasonable  because  the  “potential” 
here  is  very  smooth  near  the  lattice  points,  and  does  not 
have  a  highly  attractive  part  as  that  found  in  the  corre¬ 
sponding  solid  state’s  problem.  The  band  structure  is 
shown  in  Fig.  1  for  /= 0.375  and  r—  12.25.  We  have 
used  reduced  units,  where  kb  is  measured  in  units  of 
2 t r/c0,  and  c0  is  the  length  of  the  conventional  unit  cube 
of  the  fee  lattice.  We  have  also  repeated  the  band- 
structure  calculation  using  the  KKR  method.  We  find 
that  it  is  necessary  to  include  terms  up  to  /ma,  =  3,  espe¬ 
cially  for  large  r  values.  The  KKR  results  are  shown  by 
the  circles  in  Fig.  1,  and  are  in  excellent  agreement  with 
those  obtained  by  the  plane-wave  method.  In  addition, 
we  have  checked  that  our  results  in  the  low-frequency 
and  long-wavelength  region  are,  within  our  numerical  ac¬ 
curacy,  in  perfect  agreement  with  our  analytical  predic¬ 
tion  in  Eq.  (7). 

From  the  figure,  we  see  that  there  is  a  common  indirect 
gap  in  the  band  structure.  This  gap  is  determined  by  the 
top  of  the  “valence  band”  at  the  W  point  and  the  bottom 
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FIG.  3.  The  size  of  the  gap  at  the  optimal  volume-filling1  frac¬ 
tion  is  shown  as  a  function  of  the  relative  dielectric  constant. 


of  the  “conduction  band”  at  the  L  point.  According  to 
our  numerical  results,  this  is  true  also  for  other  values  of 
/  and  r,  as  long  as  such  a  gap  is  finite.  Note  that  in  the 
empty-lattice  limit,  the  lowest  levels  at  the  W  point  form 
a  quadruplet,  and  are  higher  than  those  at  the  L  point, 
which  forms  a  doublet.  Thus  in  order  to  have  a  gap, 
these  levels  must  be  splitted  and  shifted  in  such  a  way 
that  the  second  level  at  L  must  be  higher  than  the  lowest 
level  at  W.  We  find  that  the  gap  persists  down  to  an  r 
value  as  small  as  3,  in  contrast  with  the  value  of  7.8  re¬ 
ported  by  the  earlier  study.2  Clearly  at  a  fixed  value  of  r, 
this  gap  do  not  exist  if  either  /  is  too  small  or  too  large, 
because  in  either  limit,  the  photon  simply  sees  a  uniform 
dielectric  medium.  But  we  find  that  for  r  >  3,  as  /  is  in¬ 
creased  from  zero,  the  size  of  the  gap  first  increases  from 
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negative  to  positive,  and  then  it  reaches  a  maximum. 
Further  increases  in  /  then  decreases  the  gap  size  and 
eventually  becomes  negative  again.  Using  the  plane-wave 
method,  the  optimal  value  of  /,  /op,  is  found,  and  the  re¬ 
sult  is  shown  as  a  function  of  r  in  Fig.  2.  At  r  =  3,  /  is 
25%  and  it  decreases  with  increasing  r.  For  large  r,  we 
find  that  /op,  decreases  as  r  ~ '.  For  each  r,  the  size  of  the 
gap  at  the  optimal  volume-filling  fraction  is  calculated, 
and  is  found  to  increase  monotonically  with  r,  as  shown 
in  Fig.  3.  This  is  clearly  expected  since  the  larger  r  is,  the 
stronger  is  the  scattering,  and  therefore  the  larger  is  the 
gap- 

For  the  full  vector  wave,  the  degeneracy  factor  for  all 
the  levels  is  double  that  of  the  scalar  case.  Thus  a  much 
larger  value  of  r  is  required  to  obtain  a  gap,  and  the  lift¬ 
ing  of  the  degeneracies  due  to  the  polarizations  can  be¬ 
come  very  important.  Nevertheless,  we  can  expect  the 
plane-wave  method  to  be  also  effective  in  the  vector  case. 

Note  added  in  proof.  After  submission  of  this  work,  we 
became  aware  of  a  recent  publication  with  results  very 
similar  to  those  reported  here.  In  that  work,  results  were 
computed  only  with  the  plane-wave  method.  In  addition, 
the  transformation  of  the  band  equation  into  an  eigenval¬ 
ue  problem  was  not  used  and  thus  a  root  searching 
method  had  to  be  used  for  the  eigenvalues. 
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